Chaos-assisted free-space coupling and transmission spectra of high-Q deformed microcavities 
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We investigate the mechanism of efficient excitation of high-Q modes in deformed microcavities via the 
chaos-assisted free-space coupling. A quantum scattering theory is employed to study the free-space trans- 
mission properties, from which the coupling strength can be derived. The theory predicts asymmetric Fano-like 
resonances around high-Q modes in background transmission spectra, which is in good agreement with previous 
experimental results. The phase shift occurring in the dynamical tunneling across Kolmogorov-Amol'd-Moser 
tori is further studied, which plays a key role in the Fano-like resonance. This efficient free-space coupling holds 
potential advantages in simplifying experimental condition and exciting high-Q modes in higher-index-material 
microcavities. 

PACS numbers: 42.25.BS, 42.50.Md, 42.60.Da 



I. INTRODUCTION 

Over last two decades, optical whispering-gallery-mode 
(WGM) microresonators (or microcavities) [1] with high 
quality factors and small mode volumes have promised lab- 
on-chip applications ranging from fundamental physics to var- 
ious photonic devices, such as nonlinear optics [2-4], cavity 
quantum electrodynamics [5, 6], cavity optomechanics [7-9], 
optical communications [10, 11], low-threshold microlasing 
[12-15] and highly sensitive optical biosensing [16-20]. In 
these applications, traditionally light is coupled into the WGM 
microcavities by evanescent couplers, such as prisms [21], ta- 
pered fibers [22, 23] and angle-polished fibers [24], which 
have been validated to be highly efficient. In all of these 
coupling configurations, the microcavity is typically separated 
from the coupler by a distance of subwavelength because the 
evanescent field of WGMs extends over a very short range. 
The use of the evanescent couplers, however, is not suitable 
in some important applications. For example, a higher-index- 
material microcavity [4] can not be efficiently excited by the 
tapered fiber due to phase mismatching. In addition, the ex- 
ternal couplers degrade the high-Q property in the case of 
the over-coupling regime, and they are not convenient in low- 
temperature chambers. 

It has been demonstrated that WGMs in a specially de- 
signed deformed cavity can be directly excited by a free-space 
optical beam [8, 25, 26]. This direct free-space coupling is of 
importance because it is robust and requires less rigorous ex- 
perimental condition than the evanescent couplers mentioned 
above. The efficient free-space coupling originates from the 
breaking of rotational symmetry in deformed microcavities, 
which produces a highly directional emission assisted by the 
dynamical tunneling, different from the isotropic nature of a 
circular WGM cavity [25-27]. According to the time rever- 
sion, i.e., the reversibility of light path, free-space beams at 
certain positions are expected to couple into the high-Q modes 
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via chaos when they are on resonance. So far, this type of 
free-space coupling technique has been demonstrated experi- 
mentally to reach a resonant efficiency exceeding 50% [28]. 
A straightforward method to characterize the free-space cou- 
pling is to study its transmission property, e.g., transmission 
spectrum. In this paper, we investigate the dynamical tun- 
neling properties between the chaos and the regular modes in 
detail, and predict transmission spectra of the free-space cou- 
pling by employing a quantum scattering theory. It is found 
that the spectrum can behave asymmetrically as Fano-like 
lineshape [29], in good agreement with our recent experimen- 
tal observation [30]. With the theoretical model, the coupling 
strength to the regular high-Q mode can be further derived 
from the transmission spectrum. 

This paper is organized as follows. In Sec. II, we present 
the mechanism of chaos-assisted free-space coupling, and in- 
troduce a quantum scattering theory to predict a general trans- 
mission. It is found that the transmission spectrum includes 
the contribution from both the off -resonance background and 
the on -resonance modulation. In Sec. Ill, we study the off- 
resonance background transmission in the absence of the 
high-Q regular mode, from which the unperturbed scattering 
is characterized. In Sec. IV, the on-resonance transmission 
spectra are studied in detail. They depend strongly on (i) the 
additional phase when light travels in chaotic trajectories and 
(ii) the rate of dynamical tunneling. Sec. V rigourously ex- 
plains the chaotic states and the coupling strength, with which 
we deduce the condition of highest excitation efficiency. Sec. 
VI further investigates the phase shift occurring in the dynam- 
ical tunneling process, which plays a key role in obtaining 
Fano-resonance lineshapes in the transmission spectra. 



II. CHAOS-ASSISTED FREE-SPACE COUPLING 

Before studying the transmission spectrum, we briefly ad- 
dress the mechanism of the chaos-assisted free-space cou- 
pling. Poincare surface of section (PSOS) provides a simple 
and intuitive way to model the ray dynamics in deformed mi- 
crocavities, similar to billiards in quantum chaos. Except for 
an integrable ellipse-shaped cavity, the deformed microcav- 
ity has a mixed phase space including three types of struc- 
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FIG. 1 . (Color online) (a) False color illustration of a resonant mode distribution obtained with wave simulation. The two blue arrows denote the 
unidirectional emission toward 180° far-field direction. The strength outside the cavity is magnified for clear show, (b) Schematic illustration 
of a deformed microcavity excited by a free-space laser beam. The polar angle and the incident angle of rays x gi^s the coordinates of 
Poincare surface of section, (c) Husimi projection of the excitation state (blue) and resonant state (green) in log scale, corresponding to (a) and 
(b), respectively. The red dashed curve represents the critical line determined by sin^ = 1/n. The red dotted curves above sinx ~ 0.8 stand 
for some typical KAM tori which divides the phase space into discontinuing regions, (d) Enlarged view of the area of the dashed rectangle 
in (c). Red solid line denotes a regular mode, and black dots indicates the chaotic sea. Red arrows represent the coupling to the chaos with 
strength k, and green arrows stand for the decay into other chaotic states which contributes to the loss 71 in Eq.(4). 



tures: Kolmogorov-Arnol'd-Moser (KAM) tori, islands, and 
chaotic sea, corresponding to quasi-periodicity, periodicity, 
and chaotic motion of ray trajectories. High-Q modes are usu- 
ally localized in the regular regions [31, 32], e.g., the KAM 
tori which can be described by a perturbed Hamiltonian [33]. 
For such a localized high-Q mode, the excitation by free-space 
beam is primarily attributed to two channels: (i) angular mo- 
mentum tunneling and (ii) dynamical tunneling via chaos [34]. 
It has been demonstrated that the dynamical tunneling domi- 
nates, since the lifetime of photons that refract into the de- 
formed cavity greatly increases along chaotic trajectories [25]. 

Without loss of generality, here we consider a deformed mi- 
crodisk made of silica with refractive index n = 1.45, which 
has the boundary defined in polar coordinates as 

{i?o(l + V cLi COS* (/)), for COS0 > 
, (1) 
Rq{1 + 7] y bi cos' 0) , for cos < 
j=2,3 

where i?o and 77 represent size and deformation parame- 
ters, respectively. Cavity shape parameters are set as 02 = 
-0.1329, 03 = 0.0948, 62 = -0.0642, 63 = -0.0224. 
When 7/ = 1, a highly unidirectional far-field universal pat- 
tern of high-Q modes has been predicted [35] and demon- 
strated experimentally [36]. This emission characteristic is 
clearer by plotting the near-field pattern, as shown in Fig. 
1(a). It can be seen that there exist two major emission po- 
sitions at = 7r/2, 3tt/2, corresponding to refractive escape 
from counter-clockwise (CCW) and clockwise (CW) modes, 
respectively. Thus, we expect with a time reversed way, an 



excitation beam focused on the primary emission position at 
(f) — tt/2 shown in Fig. 1(b), can eventually excite the CW 
resonant modes. To quantitatively study this chaos-assisted 
process, we now use a quantum scattering theory to model the 
transmission spectrum, from which the coupling characteristic 
of the high-Q mode can be obtained. 

In the system consisting of a microcavity and unbounded 
medium outside, the state \tpi^) which describes the electro- 
magnetic field excited by the incident beam satisfies stationary 
Schrodinger equation 

H IV'.) = IV'.) , (2) 

where H stands for the system Hamiltonian. As mentioned 
above, not only chaos but the regular mode can also be excited 
by an appropriate free-space beam thanks to the dynamical 
tunneling. This can be demonstrated by plotting the Husimi 
projection [37] of the excitation state | ?/'.), as shown in blue in 
Fig. 1(c). Thus IV'.) can be expanded as a linear combination 
of chaos |C.) and regularity |WGM) [38], with the form 

IV'.) = a.|WGM) + J dc^'6.(w')|C.'). (3) 

The system Hamiltonian satisfies 

(WGM|ff|WGM) =ujQ- ij/2, (4a) 
(C.'|H|C.) -c^(5(cj'-c^), (4b) 
(C.|i/|WGM) =1/.. (4c) 

Here ujq and lj are the frequencies of the resonant regular 
mode and the incident light, respectively. The coupling coef- 
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ficient between jC^^,) and |WGM), governed by the dynam- 
ical tunneling, is described by V^. The decay rate 7 con- 
sists of the intrinsic loss and the chaos-assisted tunneling loss. 
In detail, the intrinsic decay rate 70 is attributed to radiation, 
material absorption and scattering losses in the cavity, while 
the chaos-assisted decay rate 71 describes tunneling into the 
chaotic states other than jC^^;). The rigorous physical mean- 
ings of the chaotic states and the decay rate will be further 
discussed in Sec. V. 

In this paper, we consider the chaotic states as contin- 
uum and use a standard quantum scattering model to interpret 
the transmission lineshape. Here we assume that \Cuj) and 
|WGM) are orthogonal [39]. Substituting Eq.(3) into Eq.(2), 
the coefficients a and b are determined by 



(wo - n/2)a^ + J duj'V*,buj{uj') = wa^, (5a) 

Vuj'a^+uj'b^iuj') ^ubujiuj'). (5b) 

On the one hand, applying a standard treatment [29], the 
coefficient b yields 



where 



w — w 



UJ — ujQ + — F{uj) 



(6) 



(7) 



The shift of resonant frequency resulted from the coupling 
is expressed as F{uj) =v.p.J k/{2tt{uj — a;'))da;', where 
v.p. denotes Cauchy's principle value. The reduced coupling 
strength n between \G^) and |WGM) is obtained through 
the Fermi's golden rule under the first Markov approximation 
[40], with 



(8) 



For high-Q mode in slightly deformed cavity whose intrinsic 
line width 7 is orders of magnitude smaller than the resonant 
frequency wo, the bounds of the integral of can be ex- 
tended to infinity, resulting in 



+00 



F{io) = ^v.p. / ' dc^'^^ = 0. (9) 

27r uj-uj' 

On the other hand, the normalization condition (?/;^/ \il>J) = 
5{ijj' — u) determines the value of a by 



(10) 

By integrating this equation over w, we have 



as K + 7 = 7t. The decay rate 74 can be measured directly 
from the taper transmission spectra. It should be noted that 74 
remains unchanged when the free-space coupling efficiency 
changes [26], which is distinguished with the fiber taper cou- 
pling. 

Finally, the transmission spectrum is calculated as 



TH = |(^J5|in)r 



|a„P|(WGM + v.p. / Auj 



(12) 



C„|5|in)|^ 



where S" is a suitable transmission operator connecting |in) 
and \G^), and | (C^j |S'|in) p describes the probability of trans- 
mitted signal [29]. To get a more general expression, we in- 
troduce a dimensionless frequency detuning defined by e = 
(w — ujo)/{k/2) and the ratio K = / k = (74 — There- 
fore, the transmission is simplified to 



Tic.) 



iK\ 



(1 + KY + £2 



|(CJ5|in)|^ 



(13) 



Here represents the crucial lineshape parameter of the 
transmission spectrum T{lj), taking the form 



TrV*{C^\S\m)' 



(14) 



where \ip^) = |WGM) + v.p. / duj "''-^'y^;' . To give a clear 
understanding, we consider two extreme cases. 

(i) In classical mechanics where the dynamical tunneling 
is forbidden, the regular mode cannot be excited. Thus there 
is no interaction between regular mode and the chaotic states 
(k 0), and the coefficients £,K oc 1/k as well as cx 
I/^/k. In this case the transmission spectra yields to 



ToH = |(C„|5|in)|, 



(15) 



which can be regarded as the unperturbed scattering. In Sec. 
Ill we will discuss the unperturbed scattering, which is of 
much concerning about the lineshape near resonance. 

(ii) When the intrinsic loss of regular mode is negligible and 
the regular mode is fully excited by a phase conjugation wave 
of its emission pattern, i.e., 7 ^ k and K ^ (or namely, 
'completely excited', which will be discussed specifically in 
Sec. V), the transmission yields a standard Fano resonance 



-|(C.|^|in)|^ 



III. OFF-RESONANCE TRANSMISSION 



(16) 



27r (^_c.o)2 + (^)2- 



(11) 



In Eq.(ll), la^^p describes the excitation probability by the 
free-space beam, from which we can deduce that the FWHM 
(full width at half maximum) of the regular mode is expressed 



Before detailing the transmission spectrum around a high- 
Q regular mode, we investigate the background scattering in 
the absence of the regular mode. It has been reported that non- 
resonant pumping in deformed microcavity can be well mod- 
eled by ray dynamics [41, 42]. In our case, the unperturbed 



4 



CD 



CD 



X3 
CD 

1 
(/) 

CD 







Radial 




Vertical 



Frequency 



2 
Displacement (^m) 



FIG. 2. (Color online) (a) Simulated transmission spectra via bound- 
ary element method, (b) Transmission spectra measured experimen- 
tally (grey dots) and fitting background oscillation corresponding to 
Eq.(20) (red solid line). These two graphs are plotted in different 
scales. 



FIG. 3. (Color online) Experimental ratio of the amplitudes of the 
transmitting and refracting channels t/r depending on tuning the dis- 
placement in radial (red round spots) and vertical (blue square spots) 
direction. Inset: SEM image of the deformed microtoroid in oblique 
view. The white arrow represents the incident beam. 



transmission is approximately determined by the interference 
between two components, according to the schema shown in 
Fig. 1(b): (i) the direct transmitted amplitude t, and (ii) the 
dissipated amplitude r via diffusing inside the cavity. To give 
a clear picture of the interference, we apply transmission ma- 
trix, and the amplitudes on the interface are indicated by 



t 9' 
9 t' 



Ein 



Eg' and Eg are related by 

Egi = a Eg, 



(17) 



(18) 



where a is a coefficient including the loss and the phase 
change in a round trip. The transition matrix element is given 
by 



S* in = — = < + 

Ein 1 — Ctt' 

= t + re'^. 



(19) 



Here r and 9 can be understood as the equivalent amplitude 
and phase shift of forward-emitted field from the cavity, re- 
spectively. Hence, the unperturbed transmission takes the 
form 

To = |(C„|^|in)p = 7-2 (l + {if + 2- cose) . (20) 
\ r r J 

In a wide frequency width 6 can be simplified as nkLoff 
with Loff representing the equivalent chaotic path length of 
the light inside the cavity. Thus the transmission spectrum 



shows periodic modulations in good agreement with numeri- 
cally simulated transmission, as shown in Fig. 2(a). In exper- 
iment, we focus the incident beam on the periphery of a de- 
formed microtoroid with the principle radius 45/im, and the 
waist of the beam is about 3^m [36]. Fig. 2(b) reveals that 
the experimentally detected transmission also oscillates peri- 
odically. The fitted ratio t/r for experiments increases with 
moving the focused spot away from the cavity in radial direc- 
tion or lifting the beam in the vertical direction, and the latter 
response is more sensitive, as shown in Fig. 3. It is because the 
silicon substrate of the device blocks a portion of the incident 
light, asymmetrically in vertical direction. Thus, the focused 
spot changing vertically results in a remarkable impact on the 
direct transmitted amplitude. 



IV. ON-RESONANCE TRANSMISSION 

We now turn to the on-resonance transmission. It is noted 
that the direct excitation probability of high-Q regular modes 
via evanescent field is negligible due to angular momentum 
mismatch, so that the amplitude of (WGM|5|in) in Eq.(14) 
has minor contribution of the transmission. In this case, the 
lineshape parameter q can be reduced to a simplified form. 
Substituting Eq.(19) into Eq.(14), we obtain the expression of 
the lineshape parameter 



v.p./du;-^K:,(a,|^|in) 
T:V*{C^\S\in) 



t/r + e*' 



(21) 



where 6 = nkLcff as mentioned above. Substituting Eq.(21) 
into Eq.(13), the on-resonance transmission can be deduced. 
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FIG. 4. (a)-(h) Calculated transmission spectra with h/2tt, ■y/2tt — 
0.003 GHz and t/r — 0.2. The phase shift between the two ampli- 
tudes varies from to 7n/4. 



In the following, we will show that the lineshape of the trans- 
mission spectrum is determined by q^, which primarily de- 
pends on 6, while the modulation depth relies on the relative 
coupling strength described by K. 

Figures 4(a)-4(h) plot calculated transmission spectra 
against 6, which experience symmetric Lorentz absorp- 
tion dips, asymmetric Fano-like lineshapes and symmetric 
electromagnetically-induced-transparency (EIT)-like peaks as 
6 varies from to 2tt. From Fig. 1(b), the transmission is 
a result of interference between two components: the direct 
transmitted light and the emitted light from the cavity, and 9 
actually describes the phase difference between them. Inter- 
estingly, the on-resonance transmission appears a symmetric 
dip on the background where the two components construc- 
tively interfere (6' = in Fig. 4(a)), while it switches to an 
EIT-like peak when they destructively interfere (6* = tt in Fig. 
4(e)). This is because when on resonance, the chaos refrac- 
tively excited by the incident beam can couple to the regular 
mode via dynamical tunneling, which results in a phase shift 
as energy couples back to the chaos. This phase shift in the 
round trip is semiquantitatively demonstrated tt in Sec. VI. 
Hence, although the background components constructively 
(destructively) interfere, such counteraction adds a destruc- 
tive modulation to the chaos, reflecting a dip (peak) on the 
transmission. 

As discussed above, the Fano-resonance transmission spec- 
tra can be regarded as the modulation of the high-Q mode 
on the off-resonance background. Such modulation depends 
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FIG. 5. (Color online) Calculated transmission spectra with 7t /27r = 
0.006 GHz and t/r = 0.2, and with K varying from to 60, for (a) 
EIT-like lineshape and (b) Lorentz-like lineshape. The background 
is set to the same level for each graph. 



strongly on the coupling strength k between the chaos and 
the regular mode according to Eq.(13). Here we study the two 
special cases: EIT-like lineshapes and Lorentz dips. As shown 
in the solid curve in Fig. 5(a), the modulation of the regular 
mode to the transmission spectrum is minor when K ~ 60. 
In this case, the excitation probability is extremely low. As 
K decreases (i.e., the dynamical tunneling is enhanced), the 
height of the FIT peak increases mono tonic ally, where the 
off -resonance backgrounds are lifted to the same level. When 
the loss described by 70 + 71 is negligible compared with 
the coupling strength k, i.e., K — > 0, the FIT peak reaches 
its maximum. Similarly, Fig. 5(b) shows that the dynamical- 
tunneling-induced dips become more obvious by enhancing 
the tunneling, as expected. 



V. PHYSICAL MEANING OF THE CHAOTIC STATE AND 
'COMPLETELY EXCITING' 

At the beginning of Sec. II, we have presented a chaotic 
state |Ccj). In this section, we will further investigate the 
meaning of the chaotic state and its coupling strength to the 
regular mode. To study the case in a general way, we expand- 
ing the chaotic state as a linear combination of an orthogonal 
set at a certain frequency, which can be regarded as uncorre- 
lated beams with each exciting chaos itself (or beam state for 
short). Using |C„(„)) to represent the normalized nth beam 
state, we have 



I a 



1 



Lj(n) 



(22) 



where «(„) stands for the corresponding weight. The coupling 
mode theory indicates 



(23) 



6 



Here ^„ and represent the electric field of the rin beam state 
and the regular mode, respectively, while g„ is the coupling 
strength between them. Thus, the equivalent coupling strength 
I Vu, I can be derived as 



IK. 



1 



'(n)5(j^ 



(24) 



As mentioned above, the total coupling strength k = 2tt\Vu:\'^. 
Considering the decaying process from the regular state, we 
have 



(25) 



where denotes the field of the chaotic state resulted from 
the regular mode spontaneously decaying via chaos-assisted 
tunneling (or namely, decaying state I'ju:)) and is the decay 
rate. Equation (23) also indicates that the decay state satisfies 



17.) = 



1 



(26) 



from which the decay rate can be deduced, as 



(27) 




0.01 



Polar angle ^(n) 



Similarly, the first Markov's approximation leads to the ex- 
pression of total chaos-assisted tunneling rate of the regular 
mode that 7c AT = 27r|7(^|^. According to Cauchy inequal- 
ity, however the regular mode is excited, the coupling strength 
K < 7CAT ■ Hence, to maximize the coupling strength, the co- 
efficients «(■„•) satisfies 



"(1) 
5(1) 



"(2) 
5(2) 



9in) 



(28) 



It is found that in this case \Cu:) = |7w), which means the 
incident light completely matches the time-reversal pattern of 
the regular mode's emission, in agreement with the case of 
'completely excited' discussed in Sec. II. 



VI. PHASE SHIFT IN DYNAMICAL TUNNELING 
THROUGH KAM BARRIERS 

Finally, we discuss the phase shift occurring in the dynam- 
ical tunneling through KAM curves [43, 44]. This phase shift 
is crucial to give rise to the Fano resonance. As shown in Fig. 
1(c), KAM tori separate the phase space into disconnected re- 
gions, between which transport is forbidden classically, but 
permitted in quantum mechanics [45, 46], known as the dy- 
namical tunneling. It is found that chaos contributes to the 
enhanced coupling strength. As shown in Figs. 6(a)-6(c), 
the Husimi projections of the excitation states inside the cav- 
ities depend on the deformation of the cavities. For ?/ — 0, 



FIG. 6. (Color online) (a)-(c) Husimi projections of the excitation 
state inside the cavity at non-resonant frequency with deformation 
parameter rj of the cavity set as [0, 0.5, 1] from top to bottom. The 
red dashed line denotes the critical line. 



the KAM tori which fill the PSOS greatly prohibit the exci- 
tation of high-Q modes. As the deformation parameter in- 
creases, regular high-Q modes are more likely to be excited. 
To evaluate the tunneling effect, we study the case that a pho- 
ton crosses an additional barrier to a regular periodic. For the 
sake of analytical expressions but without loss of the physics, 
we investigate the KAM barriers in the circular microcavity. 
The wave function of a WGM with angular momentum num- 
ber m in circular cavities takes the form 



where /m (r) satisfies the radial wave equation 



(29) 



[-V? - + 4^f^^r) . 4ur). (30) 

c r c 

Based on the stationary Schrodinger equation 

{-^V^ + V)U{r)^EU{r) (31) 

and substituted with E = huj, we deduce the effective poten- 
tial corresponding to the angular momentum number m 



V 



(n(r)2 - 



(32) 



7 



Extending m to the range of positive real numbers, from the 
classical relation m = 7r/(7r — 2%), the potential takes the 
form 



V{smx) 



2^ (tt - 2xYrl 



(33) 



where rg is the radius of the cavity. The effective potential 
V aggrandizes rapidly as x increases, and tends to infinity as 
X^ 'n 12. The tunneling coefficient across this barrier is 



1 



1 



(34) 



For the high-Q regular modes lying above sinx > 0.95, the 
effective potential V is extremely large, so that arg(ip) ^ 
—i: 12. Hence, a tt phase shift emerges in the counteractive 
process — > |WGM) |C„), which leads to the trans- 
parency of chaos, as discussed in Sec. IV. 



to the enhanced excitation of regular modes via chaos-assisted 
dynamical tunneling. A quantum scattering theory is em- 
ployed to describe the picture and to obtain the free-space 
transmission. Unlike evanescent coupling with a waveguide 
where the transmission spectra behave symmetrically, this 
model predicts three types of transmission, i.e., asymmetric 
Fano-like, symmetric EIT-like and Lorentz dip lineshapes, de- 
pending on the phase difference related to the fluctuation of 
background transmission. It is found that the Fano resonance 
is attributed to the phase shift occurring in the dynamical tun- 
neling into classical-forbidden regions. Our results provide 
a general method to evaluate the coupling strength between 
the chaos and the regular mode from the transmission spec- 
tra, which can be further extended to the quantitative study of 
the dynamical tunneling process. The efficient chaos-assisted 
free-space coupling is of importance to simplifying experi- 
mental condition and exciting high-Q modes in higher-index- 
material microcavities. 
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